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^ ■ Abstract 

, The reaction ir^d ^ nwy is calculated in chiral perturbation theory so as to facilitate an ex- 

O 
D 

Q 

QQ ' This includes loop effects in the elementary vr p — > amplitude and two-body diagrams, both 



traction of the neutron-neutron scattering length (flrm)- We include all diagrams up to 0((5 ). 



pv^j . of which were ignored in previous calculations. We find that the chiral expansion for the ratio 
> , 

Q> , of the quasi-free (QF) to final-state-interaction (FSI) peaks in the final-state neutron spectrum 

' converges well. Our third-order calculation of the full spectrum is already accurate to better than 

■ 5%. Extracting a.„„ from the shape of the entire tt d — > nn'y spectrum using our calculation in its 



present stage would thus be possible at the ±0.8 fm level. A fit to the FSI peak only would allow 
an extraction of ann with a theoretical uncertainty of ±0.2 fm. The effects that contribute to these 
error bars are investigated. The uncertainty in the nn rescattering wave function dominates. This 
suggests that the quoted theoretical error of ±0.3 fm for the most recent 7r~(i — > nn7 measurement 



^ \ may be optimistic. The possibility of constraining the nn rescattering wave function used in our 

d ' calculation more tightly — and thus reducing the error — is briefly discussed. 
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I. INTRODUCTION 



In QCD, charge symmetry (CS) is a symmetry of the Lagrangian under the exchange 
of the up and down quarks Q]. This symmetry has many consequences at the hadronic 
level, where it translates into, e.g., the invariance of the strong nuclear force under the 
exchange of protons and neutrons. However, CS is broken by the different masses of the 
up and down quarks and thus the strong interaction manifests charge symmetry breaking 
(CSB). The different electromagnetic properties of the up and down quarks also contribute 
to CSB. An important consequence of the first CSB effect (strong CSB) is that the neutron 
is heavier than the proton, since if CSB was only an electromagnetic effect the proton would 
be heavier and prone to decay. This would make our world very different, since big-bang 
nucleosynthesis is dependent on the relative proton and neutron abundances. 

While there are a number of pieces of experimental evidences for CSB 111 — including 

Q n 

recent results in dd an at lUCF m and np — > dir^ at TRIUMF j3| — one of the most 
fundamental to nuclear physics is the difference between the neutron-neutron (a„„) and 
proton-proton (a^p) scattering lengths. The scattering lengths parameterize the nucleon- 
nucleon interaction at low relative energies through the effective range expansion. In this 
low-energy region the (^5*0) phase shift 60 can be expressed as 

pcotSo = -^ + Kqp'^, (1) 

where p is the center-of-momentum (cm.) relative nucleon momentum and tq the effective 
range. This expansion is reliable for p ^ 150 MeV/c. 

Since there are no free neutron targets it is very difficult to make a direct measurement of 
the ^Sq neutron- neutron scattering length, though proposals to do so have been made. The 
latest of these suggests using the pulsed nuclear reactor YAGUAR in Snezhinsk, Russia . 
However, the more common — and so far more successful — approach is to rely on suitable 
reactions involving two free neutrons and corresponding theoretical calculations to extract 
ann from indirect data. By choosing the kinematics carefully one can detect the neutrons 
in a low-energy relative S'-wave which can be accurately described by the effective range 
expansion The currently accepted value 

a„„ = -18.5 ± 0.3 fm (2) 

is deduced from the break-up reaction nd — > nnp and the pion radiative capture process 



7r~d ^ nn'f. The first reaction is, however, complicated by the possible presence of three- 
body forces, but even after they are taken into account there are significant disagreements 
between values extracted by the two techniques. A recent nd break-up experiment reports 
Onn = —16.1 ±0.4 fm P|, i.e., more than five standard deviations from the standard value 0. 

The result of Ref. 5| is also in disagreement with another nd experiment that claims — 18.7± 

fl □ 

0.6 fm [6|. Earlier data had an even larger spread, see Ref. [3| for a review. Since the 

proton-proton scattering length is app = —17.3 ±0.4 fm (after corrections of electromagnetic 
effects |46|), there is even uncertainty about the sign of the difference The more 

negative a„„ is favored by nuclear structure calculations, where the small (but important) 
CSB piece of the AV18 potential is fitted to reproduce with ttnn = —18.5 fm |8|. 

The binding-energy difference between and ^He, which has a small contribution from 
CSB effects, is then very accurately reproduced. This would not occur were app — a„„ to 
take the opposite sign 9|. 

Because of these issues the accepted value is weighted toward the a„,„ = —18.50 ± 
0.05(stat.) ± 0.44(syst.) ± 0.30(theory) fm reported by the most recent n^d nw-f ex- 
periment Q|. The extraction in this is case done by fitting the shape of the neutron 
time-of-fiight spectrum using the model of Gibbs, Gibson, and Stephenson (GGS) ^|. This 
model was developed in the mid- 70s and explored many of the relevant mechanisms and 
the dependence on various choices of wave functions. Gibbs, Gibson, and Stephenson cal- 
culated the single-nucleon radiative pion capture tree-level amplitude to order p/M, and 
consequently ignored the pion loops that would enter at the next chiral order. Two-body 
diagrams were not fully implemented in this model. The theoretical error was dominated 
by uncertainties in the scattering wave function. Similar results were obtained in earlier 
7T~d —>■ nwj experiments carried out at the Paul Scherrer Institut (PSI) (then Swiss In- 
stitute of Nuclear Research) In the PSI experiments, only the FSI peak was fitted, 
while at LAMPF the entire spectrum was fitted. The theoretical work for the PSI results 
compared the GGS model with work done by de Teramond and collaborators jl^. The 
latter used a dispersion relation approach for the final state interaction, with a theoretical 
error of the order 0.3 fm, i.e., similar to GGS. 

In this paper we recalculate the ir^d nwj reaction using chiral perturbation theory 
(xPT). The one-body and two-body mechanisms are thus consistent and the constraints of 
chiral symmetry are respected, which is of crucial importance in this threshold regime. At 
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third order, 0{Q^), all the amplitudes of the previous calculation are included, as well as 
pion loops and three pion-rescattering diagrams. Additional advantages of the chiral power 
counting are that it gives a clearly defined procedure to estimate the theoretical error and 
provides a systematic and consistent way to improve the calculation if needed. In this first 
paper we establish the machinery necessary for a precise extraction of the nn scattering 
length. We isolate the sources of the largest remaining errors and suggest means for their 
reduction, which should make it possible to reach the desired high precision in future work. 

The paper is organized as follows. In Sec. |n] we will develop the main ingredients of our 
calculation: the Lagrangian, the explicit forms of the one- and two-body amplitudes, and a 
description of our wave functions. The numerical results are presented in Sec. IIIII together 
with our estimate of the theoretical error. We conclude in Sec. IIVI 



II. LAYOUT OF CALCULATION 

The LAMPF experiment used stopped pions, captured into atomic orbitals around 
the deuteron. The subsequent radiative decay occurs for pionic s-wave orbitals only. Thus 
the cm. and laboratory frames coincide and the pion momentum is vanishingly small. The 
neutron time-of-flight distribution of the cm. 7r~d decay width (with the photon and one 
neutron detected) can be expressed as 

dhdOs 2{2'KfM^di M^a -Ei-pi cos ^3 3 ' ' ' ^ ^ 

where ti, pi, and Ei are the time-of-fiight, momentum, and energy of the detected neutron, 
^3 is the supplement of the angle between this neutron and the photon, and k is the photon 
momentum. Here the sum is over deuteron and photon polarizations and M^^d is the mass of 
the deuteron-pion bound system, which to a very good approximation is given by the sum 
of the pion (/x) and deuteron (M^) masses. 

The matrix element Ai is the sum of four interfering parts, the quasi-free (QF) one-body, 
the one-body with final state interaction (FSI), and the two-body contributions, with and 
without FSI. These can be symbolized by the generic diagrams shown in Fig. ^ In this first 
calculation we restrict FSI to S-waves only, and subtract the plane-wave from the scattering 
wave function and include it in the QF contribution. 

We will derive the matrix elements for 'jnn — tt^c? rather than 7r~d nwj in order to 
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(a) 





FIG. 1: Generic diagrams of the contributions to TT~d nn^. Diagram (a) and (b) are the 
one-body contributions, without and with FSI. Diagram (c) and (d) are the corresponding for the 
two-body currents. The amphtudes and are described in the text. 

reduce the possibility of relative phase errors when using the 7n — > n'p amplitudes. The 
7r~d ^ nw-f decay rate can of course then be obtained by detailed balance. An explicit 
expression for the matrix element in terms of vrrf atomic ($7rd) and deuteron {(pa) wave 
functions and the pion-photon amplitude A is given by 



M.{'~)nn TT d) 



-.KM^m-^;p")A^-^{p,p), (4) 



(27r)9 VE^mK 

where p, p', and p" are the initial, intermediate, and final relative momenta of the two 
nucleons (for 77272 —>■ 7T~d), while q is the pion cm. momentum and the energy of the 
indicated particle. The meaning of these kinematic variables can also be inferred from Fig. El 
Here \E'_k(p', p) is the nn scattering wave function at total momentum — k, normalized such 
that 



(27r)3 

Its Fourier transform is given by 



klP ,Pj 



d^r 
(2^ 



- tp r 



kir,pj 



(5) 



(6) 




-k/2+p 



-k/2-p 



FIG. 2: Basic diagram for tt d ^ mi'j, defining kinematic variables. 



where \l/_k(p,r) = Z]/(2/ + l)i'^^^Pi(p ■ f), which could be compared to the plane wave 
expansion e^^ '^ = I]/(2/ + l)i^ji{pr)Pi{p ■ r). The form of vi{r) will be discussed in Sec. Ill El 
The pion-photon amplitude A can be separated into the one- and two-body amplitudes 
Ai and ^n; 



A = (27rf|:5(^)(p±^Mi(k,q)+^n(±P,k,q) 



(7) 



where the upper(lower) sign is for interaction on nucleon 1(2) of the deuteron and p = p'— p". 
The one- and two-body matrix elements are then 



A4i{'-fnn ^ n d) 



7Wn(7nn ^ n d) = J 



d^qd^p' E^/2E~i k-q , 

*^d(0; q)</'d(-q; p ± ^^)^i(k, q)*-k(p , p), 



(27r)6 ^EiE22E^ 
d^d^p'dY M^/2E~i 



(27r)9 V^i^22K 
X ^^_k(p',p). 



<l>;,(0;q)^^(-q;p"M„(±p,k,q) 



(8) 



In configuration space the matrix elements are given by 



M^, /^^i^[^^^(^) + ^^(P)^(^)]g±f(k-q).r^^^k,q)5o(r)xo, (9) 



dVt 

da 

^71 



-Mff = T.'^UO)J^ I dr^[S,u{r) + S,{v)w{r)]Af,{vXc^)v,{T)xo. (H) 
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Ml^ = E^UO)^^ I rrfr-^[5on(r) + 52(f)i/;(r)]^±(r,k,q)e±^P-^X5, (12) 

where pi^2 = — |k ± p, A^{r,'k,q) = J ^|^e~*P''".4.ii(±p, k, q), and Vo{r) is the subtracted 
scattering wave function as defined in Sec. Ill El The sums are over the two nucleons. In 
these expressions $7rd(0) = {f^ndOi)^^'^ / \fTx is the pion-deuteron atomic s-orbital wave function 
evaluated at the origin with 11.,,^ the reduced vrd mass, while S'o = ' ^'^d 5*2 (f) = 

^(3cr ■ r r ■ e)^ — cr ■ e)^) are the S- and D-wave spin structures of the deuteron, e); being the 
deuteron polarization vector. The xs's are the neutron-neutron spin wave functions for spin 
S. 

In the following subsections we will derive explicit expressions for the amplitudes and 
show how the coordinate-space wave functions are obtained. 



A. Power counting 



We start from the relativistic Lagrangian 



r - rW I r(2) , r{2). 



^(1) 



N 



-N 



il ■ Dm 
e(Ko 4 



M - ^e'^^V%"7 ■ dn 



9a 
2U 



(7 ■ D^'iT'^y^, 



8M 



-cr 



2 ^ 



1 



71 U TT , 



(13) 



where D'^ = — ieQ^A'^ and D'^"'^ = 5"^ 9^ — ieQ'!^A^ are the covariant derivatives for 
the nucleon and pion (e < 0), Qn = |(1 + t'^) and = «e"*^ are the electric charge 
isospin operators (with a/b isospin indices of incoming/outgoing pion), and kq^i = ± k„ 
{Kp = 1.793 and /t„ = —1.913) represent the nucleon anomalous magnetic moments. The 
electromagnetic field tensor is given by F^^^, = dpA^i_ — d^A^ and a^^ = |[7^, 7^^] as usual. 

The Lagrangian ()13j) is organized according to the number of powers of small momenta Q 
(here e counts as one "small" momentum). The chiral order of any graph to any amplitude 
involving nucleons, together with pions and photons with energies of order fi can be assessed 
by multiplying the Q-scaling factors of the individual units of the graph by one another. 
These factors cirG cis follows: 



Each vertex from C^^'> contributes 



• Each nucleoli propagator scales like \/Q (provided that the energy flowing through 
the nucleon line is ~ //); 

• Each pion propagator scales like l/Q^; 

• Each pion loop contributes Q^; 

• Graphs in which two nucleons participate in the reaction acquire an extra factor of 

In practice tree-level relativistic graphs must be calculated, and then expanded in powers 
of p/M in order to establish contact with the usual heavy-baryon formulation of chiral 
perturbation theory. The expressions for loop graphs that we use are also computed in 
heavy-baryon xPT, and so we do not need to employ any special subtraction schemes to 
remove pieces of loop integrals which scale with positive powers of the nucleon mass 



1J,|15|. 



Note that we will employ the Coulomb gauge in all calculations. 



B. One-body amplitudes 

There are four basic one-body diagrams, shown in Fig. El the Kroll-Ruderman (KR) 
term (a), the pion pole (b), and the s- and w-channel nucleon pole terms (c) and (d). They 
can be calculated directly from a non-relativistic reduction of the relativistic Lagrangian 
or from the amplitudes of heavy-baryon xPT (HB^PT) In addition, there are also 

pion-loop corrections at 0{Q^) as shown in Fig. EJ The loop corrections, together with the 
corresponding counterterms from C^^]^ ^ have already been calculated in the cm. frame 
to 0{Q^) for radiative pion capture on a nucleon using Coulomb gauge. The low-energy 
constants (LECs) from the third-order chiral Lagrangian were fltted to experiment, yielding 
an excellent description of the near-threshold data. 

The third-order piece of the one-body amplitudes of Ref. Q] can hence be taken as is, 
without introducing any new unknown parameters, and combined with our evaluation to 
0{Q^) of the tree-level diagrams in Fig. El using the relativistic Lagrangian (fT^ 0]. In 
order to incorporate these amplitudes in the two-body system they should be evaluated at 
the relevant subthreshold kinematics, corrected for a boost to the overall rest frame, and 
corrected for off-shell effects. These issues will all be discussed below. 
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(c) ^ ^ (d) 

FIG. 3: Single-nucleon pion-photon diagrams relevant for pion photoproduction. In this and all 
other figures a solid line represent a nucleon, a dashed line a pion and a wavy line a photon. 

^ s ^ ' s ' 

Z> ^i^^ i 




' / 

^ / 

V / 

FIG. 4: Pion loops at NNLO for — > ■nN in Coulomb gauge. Note that not all time orderings 
are shown. 

The full one-body amplitude is given (in Coulomb gauge) by 

^i(7A^ ^ ttA^) = x)zcr ■ + F2(-E'^, x)cr ■ q<T • (k X e^) + F3(_E7r, a;)«cr ■ kq ■ 

+ F4(£'^,x)icr ■ qq- e^, (14) 

where the Fj are the Chew-Goldberger-Low-Nambu (CGLN) amplitudes [3| and is the 
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photon polarization vector. The isospin channels are separated as 

where a is the pion isospin index. For 771 Ti~p this implies that Fi = v^[-^/'^^ — ■*]• 
The Fj's of Ref. Q| are evaluated with the pion energy E^^ and photon-pion cosine a: = k ■ q 
in the 7r~p 772 rest frame. In our case q = 0, -Ett = /W, and x is undetermined. Thus 
only Fi survives, the other spin amplitudes being proportional to the pion momentum. In 
charged pion photo-production Fi is dominated by the KR contribution [Fig. IHfa)]. 



1. Subthreshold extrapolation 

If the 7r~p 7n process was completely free, the CGLN amplitudes should be evaluated 
at the pion threshold E^^ = fi. However, since the proton is bound in the deuteron, the 7r~p 
energy is actually less than fi, which means that we must extrapolate to the sub-threshold 
regime. To do this we need a prescription to calculate the invariant two-body energy s^-p. 
The pion and photon energy in the 7r~p 777, rest frame can then be calculated using the 
well-known relations 

e: = '--p-'^l+^'\ (16) 

u* = (17) 

The energy available to the 7i~p subsystem, s^^-p would seem to be different depending on 
whether FSI or QF kinematics are considered. Furthermore, there are two QF situations, 
i.e., the detected neutron can originate from the one-body vertex or it can be the spectator. 
In fact, the first case is overwhelmingly favored by the kinematics of the LAMPF experiment 
and is also the one closest to threshold. The second, spectator, scenario is suppressed by 
kinematics, so even though it is further from threshold and so results in a larger shift in 
s^-p, any correction resulting from this shift is small compared to other, included, effects. 

For the QF kinematics where the detected neutron originates from the one-body vertex 
the rest frame coincides, by definition, with the overall 'ynn cm. But, in the FSI region one 
has to make a choice. The invariant energy of the n^p 777, system can be established from 

s^-p = {Md + iiy + ml-2{Md + ii)e,, (18) 
10 



where eg = ^Jra"^ + p'i is the energy of the spectator nucleon. We choose to assume that the 

spectator nucleon is on-shell and that its typical momentum Ps can be estimated through 

calculating the expectation value {pD between initial and final state wave functions. Using 

the S-state of the deuteron only, the average is given by 

P _ / drjMB -p^ + 2MVssHr)U^)do{r) 
^^'^ 4 J dru{r)U!fHr) ' ^ ^ 

We then use free kinematics for the one-body amplitudes in the QF region, and the formulas 

(fTHj) and (fT^ to calculate the energy at which the one-body amplitude should be evaluated 

in the FSI peak. The one-body amplitudes are then calculated using E* according to the 

different kinematics of the QF and FSI configurations. The theoretical uncertainty due to 

this procedure is assessed in Sec. IIII C II 

2. Boost corrections 

In general the 'jn -K'p rest frame does not coincide with the overall cm., so we have 
to adjust the Fi for boost effects. The boost corrections can be calculated by replacing the 
7?7, 7r~p rest frame kinematics by the overall 'jnn cm. kinematics in the evaluation of 
the one-body amplitudes. This changes the incoming and outgoing nucleon momenta, but 
not the photon and pion momenta. The Coulomb gauge condition e° = is retained. From 
the Lagrangian one can then deduce the following boost corrections for the reduced 
amplitudes, up to order Q'^/M'^ for the ttN reactions; 

AF<°'(£J = £| Hg.^^£ + g;) (,^ + (20) 
AF,'-'(B.) = £| ^-P~^£ + ^' . (21) 

where p„ is the outgoing nucleon momentum (= — k in the 'jn rest frame, which makes 
these amplitudes vanish). As before we have assumed Coulomb gauge, q = 0, and the 
same isospin designations as in Eq. ()15|). These corrections should thus be added to the 
amplitudes of Eqs. (fT^ and (fT3j) as given in Q|, except for an overall factor of M/47^^/s, 
which is included in the phase space in our formalism. There are also terms with new 
spin-momentum structures: 

G'»'(B,) = + - 1). (22) 
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(a) (b) (c) 

FIG. 5: The third order pion rescattering diagrams relevant for -K~d — > nn^. They are shown in 
order of expected importance as explained in the text. 

G {E.) = — (/i, - /i„ + -) j , (23) 

which will also vanish in the limit Pn — ^ — k. In the case of no n- vanishing pion momentum, 
additional terms will show up. One would expect that the first term of G^~^ should give 
the largest contribution since /ip — /in + | = 5.2 is a big number. However, because of 
the particular kinematics of the present problem, p„ ^ — k and the triple scalar product 
p„ ■ (k X e^) ^ £'^sin^3 with ^3 = 0.075. Thus, ultimately this piece of is very small 
because of the kinematics. Similarly p„-e^ ^ E^^ sin 6*3. (Additionally, only one of the photon 
polarizations can contribute.) In fact it turns out that the new spin- momentum structures 
have a negligible effect on the pion-photon amplitude and the only possible relevant boost 
corrections come from the terms in Eqs. ^K)^ and (I2H)- In the actual calculations the 
subthreshold value E* and not E^^ = fi was used in evaluating these. 



C. Two-body amplitudes 



At third order there are three pion-exchange diagrams, displayed in Fig. El Of these, the 
first one (a) is expected to give the largest contribution since its propagator is coulombic, i.e., 
behaves like 1/q^, where q is the momentum of the exchange pion 20|. This is because (for 
our kinematics) the pion energy is transferred completely to the photon and vanishes against 
the pion mass in the propagator, thus putting the pion effectively on-shell. It has been 
argued in the literature that when the intermediate nucleon-nucleon state for this diagram (as 
interpreted in time-ordered perturbation theory) is Pauli-allowed, corrections due to nucleon 



recoil need to be taken into account 



2l| . However, in our case the intermediate nucleon pair 
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is in a triplet-isotriplet state imp 
recoil correction evaluated in Ref. 



ving a relative P-wave, which is Pauli-suppressed. The 



2ll | is thus small. The second graph (b) has an extra pion 



propagator which is off-shell, reducing the magnitude of this diagram. The third two-body 
amplitude (c) has two off-shell pion propagators and is hence suppressed compared to the 
other two. More importantly, this diagram is further suppressed since in Coulomb gauge it 
is proportional to the (vanishingly small) pion momentum. 

The two-body amplitudes, corresponding to the diagrams of Fig. (a-c), are 

egA {-2iE^) 



Aua{p, k, q) 
AiibiP, k, q) 



2/. 4/2 

egA 4:iE^ 
V. 4/2 



L(P + ^)^ 



IP 



r^rlai ■ (p - j 



k+q \2 
2 ) 



[p + qj 



(24) 



TfT^(72 ■ (p + V) ■ (p - q) 

(p-^)^[/i^ + (P+^)^] 



(25) 



^llc(p, k, q) 



egA 2-E^ 



2/. 4/2 



uj{E^ - qy) 



^i-(p-^; 



'^2-(p+^; 



2 ^ 



IP 



k-q '|2 

2 / 



(26) 



where ?/ = k ■ q is the pion-photon cosine in the overall cm. In configuration space (for 
q = 0) the two-body amplitudes can be expressed as 



Aia(r,k,q= 0) 
Aiib{r, k, q = 0) 



egA -2iE^ ( 
8/3 47rr 

-^^rM / dae-^^ 



k-r I -M—CL 



T2 (T2 ■ e^e 2 



k r 



(27) 



8/3 An 

(72 ■ [(1 - a)k - + l)?]e-^(^-°)k-r 



■ r ( CTi ■ [(1 — a)k + z(/i H — )f]e 



— Q)k-r 



- -(cTi ■ e^e'(^-")''-" - (72 ■ e^e-'(5-°)k-n 
r 



(28) 

AiUr, k, q = 0) = 0, (29) 

where /i^ = a{fi'^ + u^) — oP'uP'. These expressions agree with the ones derived in Ref. [2^ 
after the sign correction of Ref. 
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D. Matrix elements 



The full matrix elements for the QF amplitudes, projected on spin-0 and spin-1 final 
states are, after taking the trace over nucleon spins and isospins 

Mo = a \el ■ e^(Fi - XF2) +kelq- e^(F2 + F3) + q ■ ej^q ■ e^fA /(pa) 



3Ci , 
72 



P2 • 4 



P2 ■ f^ji^i - + P2 ■ kq ■ e^(F2 + F3) + P2 ■ qq ■ e^F^ g{p2) + (2 -> 1) 



Ml = C ■ (e^ X e„„,)Fi - 

+ ej^ ■ (k X e„„,)q ■ e^Fs + ej^ ■ (q x e„„)q ■ e^F^j f{p' 



(k X e^)F2 



3C - 
+ ;^P2 ■ 4 [p2 ■ (e^ X e„„)Fi - (p2 ■ qe„„ - q ■ e„„,p2 + P2 ■ e„„q) ■ (k x e^)F2 

+ P2 ■ (k X e„„)q ■ e^Fg + P2 ■ (q x e„„)q • €^^4 51(^2) -(2^1), 



(30) 



IM. 



nd 



where e„„ is the polarization vector of a spin-1 neutron pair 

C = v^<l>.rf(0) 
f{p) = I r dr 



u{r)jo{pr) 



V2 



w{r)j2{pr) 



r drw{r)j2{pr). 



(31) 

(32) 
(33) 



The corresponding spin-0 FSI matrix element is easily obtained by the replacement p2 — k 
and letting 



fik) 
9{k) 



dr 



^o('^), 



drw{r)j2 I Y I fol^'j- 



(34) 
(35) 



The symmetrization (2 ^ 1) is then equivalent to an overall factor of two in the spin-0 FSI 
matrix element. Similar expressions can be derived for the two-body amplitudes and for 
higher partial waves. 

E. Wave functions 



It is possible to calculate quite accurate deuteron and nucleon-nucleon scattering wave 
functions from the well-established asymptotic states. By using data extracted from the 
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Nijmegen phase-shift analysis [23[ as well as a one-pion-exchange potential we ensure that 
the behavior of the wave function at r ^ — is correct. This yields wave functions that are 



consistent with those obtained from xPT potentials at leading order [2J|. In order to be fully 
consistent with the 0{Q^), or NNLO, operators we have derived here one should of course 
include 0((5^) corrections to the NN_ pot ential, i.e., incorporate at least the "leading" chiral 
two-pion exchange (TPE) 
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12811. This will be done in future work. 



1. The deuteron wave function 

The deuteron wave function at large distances is described by the asymptotic S- and 
Z^-state wave functions: 

u(o)(r) = Ase-''', (36) 

= + :^ + (:^) --^^ (37) 

where 7 = VMB = 45.70223(9) MeV/c [B = 2.224575(9) MeV], Ag = 0.8845(8lfm-i/2 
is the asymptotic normalization, and r] = 0.0253(2) the asymptotic D/S ratio j22|. The 
(un-regulated) deuteron wave functions u{r) and w{r) can be obtained from the asymptotic 
ones and the radial Schrodinger equation by integrating in from r = 00 



u{r) = u^^\r) -M j dr'GQ{r' ,r)[Vss{r')u{r') + VsD{r')w{r% 

POO 

w{r) = w^^\r) - M dr'G2{r\r)[VDs{r')u{r') + VDD{r')w{r')]. (38) 



Here, GQ/2{T\r) is the S- and Z^-wave Green function (propagator) and Vl'l the standard 
projections of the Yukawa OPE potential: 







Vss = 


f 




r 


Vds = 


2V2f 






Vdd = 





-fir / o 

VsD = Vds = -2V2,r- ( 1 + — 



■ 2f ^ f 1 + - + j^] , (39) 

where = 0.0750(5) is the irNN coupling constant squared js^. The coupled integral 
equation (jHH|l is solved using standard numerical techniques. 

The integrated wave functions are divergent at small distances, reflecting that short- 
range physics has been ignored. Instead of trying to model this piece as done in many 
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phenomenological A^A^ potentials, we choose to regulate it by matching with a spherical well 
solution at r = Rd- This procedure is motivated by the fundamental EFT hypothesis that 
results should not be sensitive to the behavior at small r. If it is, there are some short- 
distance physics that need to be included in the calculation, i.e., that the parametrization 
is incomplete. This hypothesis can be tested by varying the cut-off Rd over some sensible 
range. A thorough discussion of the boundary between long- and short-distance physics 
along these lines can be found in the lectures by Lepage js^. 

The D-wave (w) is matched at the boundary R^ by the continuity of the logarithmic 
derivative, which determines the depth of the well. The S'-wave is matched assuming conti- 
nuity and that the deuteron wave function is normalized to unity. The matching condition 
is then 

poo poo pRd 

1 — dru^ir) — / drw'^ir) = / dru^ir), (40) 

jRa Jo Jo 

where the left-hand side is calculated numerically and the right-hand side analytically. 

In Fig. ini these wave functions are compared to each other, to the modern chiral NLO 
wave function of Epelbaum et al. j^], and to the wave function of the Nijm93 potential 2^. 
Choosing Rd to be in the range 1.5-2 fm gives wave functions that are very close to the 
high-precision wave functions on the market. Note that the chiral wave function (thick 
dashed line in Fig. ^ deviates considerably from the potential model wave function (solid 
line) between 1.5 and 2.0 fm. We take this as an indication that short-range A^A^ dynamics 
are at play even at distances as large as 2.0 fm. Thus we vary our matching point Rd (and 
Rnn below) between 1.5 and 2 fm, where the lower limit is set because Rd cannot be reduced 
much further without the interaction becoming non- Hermit ian. It is possible to use the 
asymptotic wave functions (dotted lines in Fig. ^ in the calculation if the matrix element 
has the necessary factors of r to cancel the divergences of the wave functions as r 0. 
Such an approach is similar in spirit to the pionless effective field theory [EFT(^)] and gives 
analytic expressions for the matrix element. 
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FIG. 6: The S- (top panel) and D-state (bottom panel) parts of the deuteron wave functions. Our 
wave functions are labeled by the choice of matching radius -R^. These are compared to the NLO 



chiral wave function of Ref. 



23] and Nijm93 



33] 



2. The nn scattering wave function 

The scattering wave function can be calculated in a similar way. However, the asymptotic 
state is now described by the phase shift according to 



^„k(r,p) ~ 



gf5o gin(pr + 6q 



(41) 



r pr 

where 5q is calculated from Eq. with given values of p, anni and tq. In the limit of 
vanishing momentum this wave function reduces to 1 — a/r as it should. If higher partial 
waves can be neglected the only integral equation we need is the one for the ^5*0 channel, 
which is 

v^{r) = v^^\r) - M H dr'G,{r' ,r)Vss{r')vo{r'), (42) 
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FIG. 7: Our ^5*0 nn scattering wave functions (thin lines) at p = 10 MeV/c and for varying as 
indicated. Here ann = —16 fm and rg = 2.8 fm. A comparison is made with the wave functions of 
the GGS model (thick lines). The latter are explained in Q], from which reference the figure was 
adapted. Reprinted figure with permission from W. R. Gibbs, B. F. Gibson, and G. J. Stephenson, 
Jr., Phys Rev. C 11, 90 (1975). Copyright (1975) by the American Physical Society. 



where GQ{r', r) is the free S'-wave two-body propagator. This wave function is regularized 
at r = Rnn by matching the logarithmic derivative of a spherical well solution. The depth of 
the spherical well is hence energy dependent, a treatment related to the energy- dependent 
potential used by Beane et al, j^J]. As for the deuteron we vary Rnn between 1.5 and 2 fm 
in order to test our sensitivity to short-range dynamics. 
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In Fig. [7| our nn wave functions are plotted together with the wave functions used by 
GGS. The wave functions are quite similar in that they tend to the same asymptotic limit 
(the curve labeled Zero Range) for larger r. Thus they are very close to each other for 
r ^ 1.5 fm. There are, however, a few important differences between the nn wave functions 
in the two calculations: Firstly, the GGS wave functions have been derived using the Reid 
soft-core potential (RSC) (with the old larger value for the irNN coupling constant) for the 
long-range part, while we used one-pion exchange only. This explains the slight difference in 
the size of v{r) at r = 1.4 fm. Secondly, GGS match at a fixed value of Rnn = 1-4 fm, while 
we vary Rnn- Thirdly, our wave function uses a spherical well solution (sin/tr) and GGS 
assume a polynomial of fifth order, where the magnitude and first two derivatives vanish 
at r = and are matched to the RSC solution at r = Rnn- The assumption of vanishing 
derivatives is equivalent to using a hard-core potential at short distances, while many chiral 
potentials have a softer behavior for small r (see, e.g., Ref. 27|). [The different shapes of 
the GGS wave functions were obtained by adding an extra term r/r^(r — i?„„)^/pr to the 
short-range piece of their wave function.] The combined effect of all this is that our wave 
function has most variation around r = 1 fm, while the GGS wave functions varies most 
around ~ 0.7 fm. As we will see later, these differences have a strong influence on the 
assessment of the size of the theoretical error in the extracted a„„. 

An obvious improvement of our calculation would be to use wave functions whose short- 
range behavior is constrained by other observables, thus reducing the uncertainty. We could 



also compare to wave functions of modern high-precision potentials, e.g., Refs. 
the recent N'^LO chiral potentials Another extension would be to include higher 

partial waves. 

In the actual calculations we subtract off the plane- wave 5*- wave contribution jo{pr) from 
the scattering wave function [vq = vq — rjo^pr)] and then calculate the full plane- wave (QF 
or PW) contribution using 6*^ "" without partial-wave decomposition. 
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FIG. 8: Time-of-flight distribution for the Tr~d — > nwy decay rate. This and all following spectra are 
calculated assuming a„„ = — 18 fm, tq = 2.75 fm, and ^3 = 0.075 rad, unless otherwise indicated. 
The plot shows the contributions from the LO KR (dotted line), NLO one-body (short-dashed 
line), NNLO one-body (long-dashed line), and NNLO one- and two-body (solid line) amplitudes. 
The two peaks are labeled QF and FSI from their dominant contributions. 

III. RESULTS 
A. Convergence 

The calculated differential decay width is shown in Fig. |H1 for the LO KR term only and 
with the NLO and NNLO one- and two-body amplitudes added in succession. The spectrum 
shows two separate peaks, labeled QF and FSI from the dominant contributions that give 



20 



rise to them. It is clear that the LO curve is very similar to the full calculation and that 
the corrections of higher orders mainly affect the magnitude, but do have some impact on 
the shape. The evolution is most easily assessed by forming the QF to FSI peak ratio at 
the various orders. At LO the ratio is 2.58, at NLO 2.59, at NNLO one-body 2.60, and at 
full NNLO 2.71. Thus the LO, NLO, and NNLO one-body results are very close to each 
other, but do not contain the full dynamics that the NNLO two-body amplitudes provide. 
For a good quantitative result it is important to include the full NNLO amplitude — using 
only one-body amplitudes would give a wrong answer at this order. 



B. Sensitivity to a„„ 

In Fig. El the decay rate is plotted for various choices of a„„. The curves have been 
rescaled to coincide at the QF peak, to facilitate comparison of the relative height of the 
QF and FSI peak. This is done since in the LAMPF experiment Q| the scattering length 
is extracted by fitting the shape of the spectrum, not the magnitude of the decay rate. Note 
that only the height of the FSI peak changes, the valley between the two peaks is largely 
unaffected by the value of a„„. 

The theoretical error in the extraction of the scattering length has several sources and they 
will be investigated and estimated in the following paragraphs. The error in the extracted 
CLnn can be related to the error in the decay rate F by 

AT _ dT ann^ann / 

F rl F ' V / 

where the actual calculations (Fig. IH)) give that -^^^^ = 1-21 at Unn = —18 fm. Thus 

Aa„„ AF 

— ^ = 0.83—, (44) 
a result we shall use repeatedly in what follows. 



C. Theoretical error bar 



We estimate the theoretical error under the assumption that the entire time-of-flight spec- 



trum is fitted. To the best of our know 
FSI peak, which limits the kinematics 



ed, 



11 



e, previ( 



e, prev ious work have only considered fitting the 
Because of the large relative momentum 
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FIG. 9: Time-of-flight distribution for the TT~d nwy decay rate. The spectrum is calculated for 
different choices of Cnn as indicated. The curves for a„n = —20 and -16 fm have been (slightly) 
rescaled to coincide with the QF peak of the ann = — 18 fm curve. 

in the QF region, this extended analysis will have significant importance for the size of the 
error. We will use a nominal value of ann = "IS fm in our estimate of the error. 

1. Neglected higher orders in the 7T~p jn amplitude 

The present calculation ignore pieces of the ii~p — >■ jn amplitude of 0{Q^) or higher, 
which is thus three orders down from the leading piece of 0{Q). One might think that 
the error would then be of the order (u/A)^, since the first dynamical effects not explicitly 
included in our Lagrangian are associated with A-isobar excitation and so the 'high'-energy 
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scale is A, the A-nucleon mass difference, rather than A^. This is supported by the results 
of Ref. [l8| where the fitted counter terms had unnaturally large coefficients when expressed 
in units of GeV~^. However, in Ref. Q] the 0{Q^) one-body (single-nucleon pion photo- 
production) amplitude was fitted to actual data for = 142 MeV and higher (roughly 
10 MeV above threshold). The error in our calculation is thus introduced only in our 
extrapolation of the amplitude to a subthreshold energy, denoted by u*. Compared to the 
leading 0{Q) term, this gives a correction (cJq — u;*'^)/A^ ~ 4%. This is a special, very 
beneficial feature of the pion absorption process: since the pion momentum is vanishing 
there is no angular dependence and the amplitudes depend only on the photon energy. This 
error should include the errors due to uncertainties in the LECs fitted in A simple 

calculation based on the LEG fit errors and the formulas for the CGLN amplitudes gives 
corrections of the order 3% or smaller, which is in line with the above 4%. 

Since the extrapolation photon energy is roughly the same at the QF and the FSI peak, 
this error should add roughly equally to both peaks, which reduces the error on the neutron 
time-of-flight spectrum to less than the ~ 4% estimated above, since now only the shape is 
fitted. The actual calculations confirm this: the spectrum using extrapolated amplitudes (see 
Sec. Ill BT]) differs by only 1.1% in the FSI peak from the spectrum with amplitudes evaluated 
at threshold. The corresponding error in ann is thus 0.95% or 0.17 fm for ann = —18 fm. 

2. Boost corrections 

The contribution of the boost corrections [Eqs. ()20|) and (I2U)] is of the order /i^/2M^ ~ 1%, 
but the change occurs in the same direction in both peaks. Thus the relative change is much 
reduced (0.14% in Ffsi/Fqf, i.e., 0.11% or 0.02 fm in a„„) and can be completely neglected 
for the present purposes. After rescaling, the boosted curve cannot be distinguished from 
the original one. We use the calculated 0{Q^) boost correction of 0.14% as a conservative 
estimate of the boost error introduced at higher orders. The boost correction is included in 
all plots. 
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3. Off-shellness 



In calculating the one-body amplitudes we tacitly assumed that both nucleons were on- 
shell. This introduces an "off-shellness" error that should be estimated. It is well known that 
field transformations can be employed to trade dependence of the one-body amplitude on 



the "off-shellness" of the nucleon, P ^ two-body amplitude 

follows: using field transformations such as those employed in Refs. 



371 38 1 . This is done as 
391 . |40| the dependence 

of the 7n — » 7r~p amplitude on the nucleon energy is replaced by dependence on ^ plus 
terms in C\l^ and beyond. This means that the one-body amplitude for the photoproduction 
process now has no "off-shell ambiguity", although we do acquire additional pieces of the 
two-body amplitude for the charged-pion photoproduction process. (This argument is shown 
graphically in the sequence of diagrams in Fig. ^|) The new contribution, depicted in 
Fig. Unfc), involves a 77r7r vertex from d-^li, and so is 0{Q^). This two-body effect is thus 
p^/M^ ~ /i^/M^ ~ 2% down from the NNLO two-body diagrams, which contribute 7.1% 
to the rescaled decay rate (according to Fig. |HI). Consequently the error in ann from any 
potential "off-shell ambiguity" is approximately 0.02 x 0.071 x 0.83 = 0.12% or 0.02 fm. 



4. 0{Q'^) two-body pieces of the amplitude 

A larger effect comes from 0{Q^) two-body pieces of the '-/nn ^ 7r~d amplitude, such 
as the one depicted in Fig. ^2 A naive estimate indicates they should be ~ p/A^ ~ 
yu/A^ ~ 20% of the 0{Q^) two-body diagrams. This estimate is supported by studies of 
pion photoproduction to 0{Q^) in xPT j22|. This suggests that 0{Q'^) two-body effects are 
roughly a 0.7% effect in a„„. 



5. Error from wave functions 

By changing the matching points and Rnn between 1.5 and 2.0 fm, we tested the error 
introduced by our ignorance of short- distance physics in the A^A^ wave functions. This change 
was significant for the nn scattering wave functions, as shown in Fig. El The resulting error 
in ann turns out to be —0.6 fm (3.3%) or smaller. A similar spread was obtained with wave 
functions calculated from "high-quality" A^A^ potentials, e.g. Nijm I and Nijm II 

Q- We 

note that both these potentials have x^/d.o.f. = 1.03 with respect to the 1993 Nijmegen 
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(a) 



(b) 




(c) (d) 

FIG. 10: The off-shell nucleon in (a) can be taken care of by extracting a meson exchange from the 
deuteron wave function as in (b). The off-shell part absorbs the closest propagator and becomes 
the two-body diagram (c), splitting off the on-shell amplitude (d). The cross indicates an on-shell 
nucleon. 




FIG. 11: Typical two-body operator at 0{Q^). The sliced photo-nucleon vertex is from C^^jif. 

database and identical nn scattering lengths. They differ only in their treatment of the 
heavy mesons, indicating that our calculation is sensitive to truly short-range parts of the 
A^A^ interaction. Our results with the Nijml and Nijmll potentials also suggest that this 
short-range sensitivity has a greater impact on the extracted a„„ than does our neglect of 
two-pion exchange. We are confident that this uncertainty could be considerably reduced 
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FIG. 12: Time-of-flight distribution for the 7r~(i — > nnj decay rate. The spectra are calculated 
for ann = —18 fm and different choices of Rnn as indicated. The case of Rnn = 2.0 fm and 
o-nn = —18.6 fm is also plotted. 

by finding other observables that constrain the wave function, in particular its short-range 
behavior. This will be pursued in future work. Note that the change in Rnn not only changes 
the height of the FSI peak, but also the valley region. This feature could potentially be used 
in a fit procedure to distinguish the Rnn dependence from a change in a„„. 

Indeed, if one focuses only on the FSI peak then the variation in the spectrum due to the 
use of different wave functions is significantly smaller than the one discussed in the previous 
paragraph. If we adjust both calculations to agree in the valley region we find that the FSI 
peak height only differs by 0.6%. (See Fig. UHl) This corresponds to an uncertainty in ann 
of ±0.1 fm. 

Meanwhile, effects due to the bound-state wave function chosen are also small. Changing 
the deuteron wave function by varying R^ from 2.0 fm to 1.5 fm would alter the extracted 
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FIG. 13: Time-of-flight distribution for the tt d ^ nri'j decay rate in the region of the FSI peak. 
The spectra are calculated for a„,„ = — 18 fm and different choices of Rnn as indicated. 



O'nn by 0.55% or 0.10 fm. Using the Bonn B deuteron wave function instead of the EFT- 
motivated wave function yields a Aa„„ of 0.56% or 0.10 fm. 



6. Higher partial waves 

The error from neglecting higher partial waves in the rescattering wave function can be 
estimated in the following way. The higher partial waves are only substantial for large 
relative energies and are thus negligible in the FSI peak region. In the QF peak, the relative 
nn momentum is roughly 80 MeV/c, which means that the S-wave phase shift is 6q ^ 60°, 
while the P-wave phase shifts are typically 6i ^ 5°. The P- to S-wave amplitude ratio 
can then be estimated as Ai/Aq ~ sin^i/ sin(5o = 0.10. From Fig. El the S-wave FSI 
amplitude at the QF peak is Aq = VO.OSO = 0.17 and thus the P-wave FSI amplitude is 
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FIG. 14: The decay rate separated into spin-0 and spin-1 contributions. The subtracted FSI 5*0 
partial wave is also distinguished. 

Ai ~ O.IOAq = 0.017. However, since the P-waves are spin-1 and the 5'-waves spin-0 and 
the two do not interfere, the influence of the P-wave should be related to the QF spin- 
1 amplitude, which is Bi = VO.295 = 0.543. The error in the calculated QF peak is then 
2 1 1 1 Pi 1 1 cos 6 1 , where 6 is the unknown phase angle between Ai and Pi . Using the maximal 
possible error (setting cos ^ = 1) seems overly pessimistic, so we instead choose the average 
(I cos^l) = 2/tt. The relative error at the QF peak is then ^|y4i||Pi|/FQF ~ 2.9%, yielding 
an error in the extracted scattering length of 2.4% or 0.43 fm. This should be regarded as a 
conservative estimate of the error of neglecting P-waves for two reasons. Firstly, the P-waves 
could interfere destructively with each other. And secondly, we implicitly assume that the 
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radial integral for P-waves is of the same magnitude as for S'-waves, whereas it is probably 
smaller. Most importantly, it is possible to actually calculate and include the P-waves. This 
error can thus easily be pushed to higher partial waves and so made substantially smaller. 
This will be done in future work. 

Note that this P-wave error is much larger than the one estimated by GGS The 
reason is that the P-waves only contribute at large relative energy, i.e., under the QF peak, 
where they can interfere with the QF amplitudes, thus changing the QF to FSI peak ratio 
and the extracted a„„. The GGS error estimate assumes that the nn opening angle is smaller 
than 30°, which restricts the kinematics to the FSI peak region only and thus does not apply 
to the entire range of neutron energies used in the LAMPF extraction. Also in the work of 
de Teramond et al. 3] as used by the PSI group |l2|, only the FSI peak region is fitted, 
which gives a small P-wave contribution. Thus, as far as we can ascertain, our analysis is 
the first that estimates the interference of the FSI P-waves with the spin-1 QF amplitude. 

If this effect was not included in the analysis of the data in Ref. then the Aa^n of 
approximately 0.43 fm we have found here should be included in the theoretical uncertainty 
quoted in that work. However, correspondence with one of the authors of Ref. Q| suggests 
that FSI in A^A^ P-waves was included in the version of the GGS model used for the 



extraction there 
value for a„r,. 



4l| |. This source of uncertainty would then not be present in Ref. [10|'s 



7. Sensitivity to tq 

An estimate of a change in ro due to CSB can be obtained by assuming that the relative 
change in tq is similar to the relative change in flTVAr. Thus ~ so that Arg = -^Aa = 
^jl^l.S = 0.23 fm. The sensitivity to the effective-range parameter ro was tested by varying 
it away from its nominal value 2.75 fm, using a conservative spread of ±0.25 fm. This 
changes the FSI peak by 1.4% (after rescaling to the QF peak) and thus indicates a change 
in the extracted a„„ of 1.2% or 0.21 fm. On the other hand, the error suggested by analysis 
of different experimental determinations of tq is ±0.11 fm j7|. If tq is instead varied over this 
narrower range the resultant Aa„„, is only 0.5% or 0.09 fm. We will use the latter, smaller, 
error in our error budget. 
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TABLE I: Error budget for the extraction of a„„ from the 7r~(i — > nri'y reaction as it was performed 
in Ref. The calculation of the absolute errors assumes a scattering length of -18 fm. The total 
error is summed in quadrature. 



Source 


Relative error (%) 


Absolute error 


Off-sheh 


0.07 


0.01 


Boost 


< 0.11 


< 0.02 


Subthreshold 


0.95 


0.17 


0(g4) 2B 


0.7 


0.12 




0.5 


0.09 


Dep. on Rd 


0.55 


0.10 


p-wave in FSI 


< 2.4 


< 0.43 


Dep. on Rnn 


< 3.3 


< 0.60 


total 


< 4.3 


< 0.78 



(fm) 



D. Error budget 

The errors are summarized in Table HI The first four errors are due to uncertainties in 
the amplitudes, while the last four are due to the wave functions [3] • We consider the total 
error of < 4.3% to be a very conservative estimate. 

Note that if a„„ is extracted only from data in the FSI region then the last two errors 
drop to 0.2% and 0.5% respectively, while a number of the other errors listed in Table IJ are 
also reduced. We find that an extraction performed using only data from this section of 
the neutron time-of-fiight spectrum would have a theoretical uncertainty of ±0.2 fm. This 
confirms the conclusion of GGS from thirty years ago. The significantly reduced theoretical 
uncertainty comes at a price though: one must sacrifice the large number of counts acquired 
under the QF peak. We have argued above that the last two errors quoted in Table |l] can 
be decreased by additional theoretical work on radiative pion-capture on deuterium, and 
therefore we hold out hope that in future a xEFT extraction of a„„ which has an accuracy 
of ±0.3 fm (or better) and uses the full neutron spectrum obtained in Ref. jl^ can be 
performed. 

One reason for this optimism is the convergence of the chiral expansion for this reaction. 
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which can be made more exphcit by computing the QF to FSI peak ratio for the different 
orders. From Fig. |H1 we obtain 



Tqf 



(2.580 + 0.014 + 0.112 ± 0.039)(1 ± 0.05), (45) 



Tpsi 

where the first parenthesis contains (in order) the contribution of the LO, NLO, NNLO, and 
the error in the chiral expansion. The second parenthesis shows the error due to effects in the 
wave functions. Note that modifying the wave functions by including two-pion exchanges, 
P-waves, or different short-distance dynamics would already change the LO calculation, 
which is why we choose to write this error as an overall factor. The smallness of the NLO 
and NNLO one-body terms can perhaps be an effect of the particular kinematics of the 
present problem, especially that the pion momentum is vanishing. On the other hand, the 
comparatively large NNLO two-body contribution is most likely a result of a combination 
of two effects: Firstly, the two-body currents allow for momentum sharing between the 
nucleons, which would be of importance in the QF region. Secondly, in the leading two- 
body diagram [Fig. El^a)], the coulombic propagator was power-counted as l//i^. However, 
because of the small deuteron binding energy, the typical momentum is instead of the order 
7 = ^/MB = 45.7 MeV 



42| |. Since 7 <^ yU this further enhances this diagram. 



IV. CONCLUSIONS 

In this paper we have calculated the vr^rf — ^ 77,77,7 reaction, using xPT pion-photon ampli- 
tudes and EFT-inspired wave functions. The errors in the extracted scattering length from 
the operators are of the order 1%. These errors include effects that were not considered 
by Gibbs, Gibson, and Stephenson (GGS) e.g., errors from extrapolating the single- 
nucleon amplitudes sub-threshold, the boost of the 771 — > n'p amplitude from the 771 rest 
frame to cm., the effects of off-shell nucleons, and more complicated two-body mechanisms. 
A key improvement is that we have included the full two-body amplitude at third chiral 
order and have found that on the scale of the other errors it has a substantial influence on 
the extraction of the scattering length. 

Nevertheless, if a„„ is extracted from the FSI region alone our analysis within xEFT 
confirms GGS's result for the theoretical uncertainties, putting them at ±0.2 fm. On the 
other hand, if — as was done in the most recent a„„ extraction 10| — the entire shape of the 
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neutron spectrum, including both the QF and FSI peaks, is used for the extraction, then 
the uncertainty in the scattering wave function at small distances and the neglect of higher 
partial waves is a potentially large source of errors, maybe as large as 4.3%. This might 
seem like a large uncertainty, since it is almost three times larger than the 1.5% estimated 
by GGS. But, as was argued in Sec. 1111 C 51 some of the assumptions behind their error 
estimate do not seem to apply for the entire kinematic range spanned by the data from the 
LAMPF experiment. This tempts us to suggest that the error estimate given in Ref. is 
optimistic and should be increased. 

We plan to improve our model in the near future by constraining the short-distance part 
of the nn wave function using other observables and by incorporating higher partial waves. 
We will report these results in a future publication. We also plan to fold our model with 
the neutron detector acceptance and the experimental geometry in order to extract the nn 
scattering length from the data of Ref. . 

Overall we conclude that the 7r~d ^ nn^y reaction has some very desirable features that 
makes it extremely suitable for the extraction of the neutron-neutron scattering length. The 
vanishing pion momentum obviously favors a xPT calculation and also reduces the number 
of contributing terms dramatically, leading to the dominance of the Kroll-Ruderman term. 
The fact that the extraction is done by fitting the shape of time-of-flight spectra rather than 
an absolute decay rate reduces many errors further still. 

The reaction '-yd —>■ nmi^ could be used as an alternative and complementary way to 



extract the neutron-neutron scattering 
.ev.w Q a„d >ate. pape., e.,, 



nth. This reaction has been considered before, see 
A chiral calculation should be feasible for this 
reaction, and could benefit from the work of the present paper. With a threshold photon 
laboratory energy of 149 MeV it should be accessible at existing experimental facilities, 
e.g., Hl7S@TUNL after the planned upgrade and MAX-lab in Lund, Sweden. After the 
submission of the manuscript, a calculation of 76? — >■ nnir^ using chiral perturbation theory 
along lines similar to ours has become available 
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